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ТРЕХМЕРНАЯ МЕХАНИЧЕСКАЯ МОДЕЛЬ ДЛЯ ОПИСАНИЯ 
БОЛЬШИХ УПРУГИХ ДЕФОРМАЦИИ ПРИ ОДНООСНОМ РАСТЯЖЕНИИ 
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Предложена трехмерная механическая модель для описания больших деформаций нелинейно-упругих тел. 
Модель отражает различные случаи анизотропии: ортотропный, трансверсально-изотропный, а также изо- 
тропный материалы. Рассмотрено одноосное растяжение. Изучены две характеристики нелинейного мате- 
риала: аналоги модуля Юнга и коэффициента Пуассона. Анализируется влияние внутренних параметров 
модели на характеристики материала. 

Ключевые слова: нелинейная упругость, анизотропия, математическая модель, модуль Юнга, коэффици- 
ент Пуассона, одноосное растяжение. 


Введение. В теории упругости предложены многочисленные определяющие соотношения нели- 
нейно-упругих тел для случая больших деформаций. Основные из них (Синьорини, Муни-Ривлина, 
Мурнагана) подробно разобраны в [1]. В их основе, как правило, лежат представления удельной 
потенциальной энергии деформации через инварианты тензора деформации. Ниже рассмотрен 
другой подход: моделирование соотношений связи «напряжение — деформация» с помощью ме- 
ханической конструкции. 

В работе [2] была предложена двумерная математическая модель (ММ) для описания ме- 
ханических свойств полимерных материалов. Данная ММ имеет различные естественные возмож- 
ности для развития. В настоящей статье подробно рассмотрена трехмерная ММ в упругой поста- 
новке. 

Вывод соотношений новой ММ сделан на примере одноосного растяжения. Модель пред- 
ставлена пространственной механической (стержневой, пружинной) конструкцией с разными гео- 
метрическими (длины) и механическими (коэффициенты жесткости-упругости) характеристиками 
внутренних связей. При этом рассматриваются взаимозависимые свойства продольного и попе- 
речного деформирования. 

В литературе обычно большое внимание уделяется модулю упругости Юнга как основной 
характеристике свойств материалов. Тем не менее, несмотря на малую область изменения, коэф- 
фициент Пуассона является одинаково равноправной с модулем Юнга механической характери- 
стикой [3]. Современные исследования [4, 5] выявили его связь с различными важными физико- 
механическими свойствами материалов. Определение коэффициента Пуассона напрямую из экс- 
перимента является сложной задачей. Вычисление же его по формулам линейной теории упруго- 
сти с использованием любой пары из констант материала (Еи К; ми К, Еи И) приводит к большим 
погрешностям, если не обращать внимание на необходимую большую точность исходных данных. 
Кроме того, усложняет задачу еще и то, что эксперименты проводятся на разных образцах и при 
разных условиях, что не может не сказаться на конечном результате [3]. 

Коэффициент Пуассона до сих пор не объяснен до конца с точки зрения механизма де- 
формирования. В данной работе предложено возможное направление анализа этого механизма. 

Прототипом предлагаемой модели можно считать и модель кристалла, описанную 

Р.Фейнманом [6]. Автор не ставил задачи анализа больших деформаций, а определил упругие 
постоянные кубических кристаллов через энергии взаимодействия атомов. 
Основные положения модели в упругой постановке. Рассмотрим элементарный прямо- 
угольный параллелепипед упругой анизотропной сплошной среды с заключенной в него трех- 
мерной механической конструкцией (рис.1). Ее узлы привязаны к центрам граней параллелепипе- 
да. Конструкция представляет собой систему стержней (пружин), отражающих свойства связей 
между гранями. «Пружины-связи» данной модели считаются линейно-упругими. При этом модель 
является трехмерной, и особую роль в этом случае играют диагональные связи. 
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Рис. 1. Геометрия связей модели 


Характеристики жесткости связей рассматриваются как некоторые интегральные оценки 
силовых взаимодействий общей структуры материала, но не отдельных атомов в отличие от [6]. 
При этом допускаются и отрицательные значения жесткостей связей, что является отражением 
реальных сил притяжения в структуре материала. 

Если длины связей А.А. АзА., АзАв не равны, то мы получим «геометрически» анизотроп- 
ную модель. Если не равны коэффициенты упругости разных связей, то такую анизотропию ус- 
ловно назовем «механической». 

Выбором различных значений геометрических и физических параметров можно придать 
модели сложную анизотропию. При этом количество параметров соответствует представлениям 
теории анизотропных сред. В общем случае анизотропного материала необходимо задать 21 па- 
раметр модели: 18 механических параметров — коэффициентов упругостей связей и 3 геометри- 
ческих параметра — линейные размеры модели. Это согласуется с общими положениями анизо- 
тропной теории упругости, где вводится 21 константа упругости материала. 

Уравнения модели для одноосного растяжения ортотропного тела. Рассмотрим случай 
растяжения «механически» и «геометрически» ортотропной модели. Обозначим длины граней эле- 
ментарного прямоугольного параллелепипеда (длины связей) А.А’ =2а, АА. =2, АА =26у. 


Диагональные связи между @0 и Бу обозначим /›, между № и су-Пу, между ау и су- Ру. Коэд- 
фициенты упругости связей обозначим соответственно К, К», К., К, К и К г 
Рассмотрим одноосное растяжение модели силой и. ‚ приложенной в вершинах А; и Ав 


вдоль связи @(. Тогда после деформации длины связей станут соответственно А, В, С [1 МиР 
причем связи А, [и Ррастянутся, а В, Си М- сожмутся: 

А=а+0., В=Б-д,, С=с-д., [=Ю+д, М=т-д,, РЕРо+д,, (1) 
где 0, ,5ь,0,., б/, 0 


У О — знакоположительные удлинения соответствующих связей. 


Выпишем соотношения для длин диагональных связей в отсчетной конфигурации: 


2 2 2 2 2 2 2 2 2 
Ч +5 — Г . Ч + Со = Ро . Ь + Сб — По (2) 
ИВ актуальной конфигурации 
А+В?=Ё, А’+С”=Р”, В" +С*= М. (3) 


Учитывая (2) и подставляя (1) в выражение (3), получаем: 
2а0, +56,? = 2650, + = 216) + 62 , 2406, +6, — 206, +6,?= 206, +5’, 


266, +6, —26,+6,7= -2,6, +6... (4) 
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Для деформированных углов (рис.2) получаем: 


























@0 + ба 
с05ф=— = . 
1 +9, 
ш@=—= И в 
Е 1+9, 
т (5) 
со5у =— = ) 
№ т-95, 
ши =— = 5 , 
по -9, 
соз@ =^ = №0 ыы , 
Рот 6, 
. т0=— = в, 
Рис. 2. Углы в актуальной конфигурации Ро+ 6, 


Составим уравнения статики в узлах в проекциях на три ортогональные оси (рис. 3): 
Е, =К, +2, с0озф+2К,с0$0, 
2, т +2, созу + К, =0, (6) 
—2® зшф+2А, зшу + А, =0, 
где К, - реакция соответствующей связи. 
Физический закон для каждой связи: 


В, =©д, ‚ 1=а,б,с1,п,р. (7) 











Рис. 3. Внешняя сила и реакции связей модели 


Подставив в (6) выражения (5) и (7), получим нелинейные уравнения: 














о А И Но В 
№ + 1 Ро+дь Ро + р пол 
20, ть. 2,6, 4% _ ,д, =0. (8) 
№ +5, т-б„ 


149 


Физико-математические науки 








Уравнения (4) и (8) образуют систему для определения неизвестных 
ЗО, зоо Он б,. Эта система является математической моделью, соответствующей 


предложенной трехмерной механической конструкции связей. 
Задавая деформацию 0, в качестве исходной величины, из (4) и (8) можно определить 
функции: 
Е, (0) > бь(да), 6. (ба), 51(0,), 6,(6.), 6,(0а) : 
Из зависимости Р.,(б,) можно определить нелинейную характеристику, связывающую 


растягивающее напряжение с продольной деформацией, являющуюся аналогом модуля Юнга ли- 
нейной теории упругости. Аналогично из зависимости д,(0,) получаем нелинейную характери- 


стику связи поперечной и продольной деформаций, являющуюся аналогом коэффициента Пуас- 
сона линейной теории упругости 


Задавая силу РЁ, в качестве исходной величины, можно определить зависимости 
5,(ЁЕ,), 9ь(Е,), де(Е,), 9(Е,), 9,(Е,), 6,(Е,) и далее характеристику, аналогичную 
податливости. 

Таким образом можно рассмотреть одноосное растяжение модели силой Р», приложен- 


ной в вершинах А; и А» вдоль связей Ь, и одноосное растяжение модели силой Г’. , приложен- 


ной в вершинах Аз и А. вдоль связей Су. При этом структура уравнений сохраняется, меняются 


только индексы. 
Уравнения модели для одноосного растяжения трансверсально-изотропного тела. Рас- 


смотрим трансверсально-изотропный материал. Пусть в отсчетной конфигурации су = В ‚ тогда 


Ро =. Найдем две характеристики материала: «коэффициент Пуассона» и «модуль упругости 


Юнга» в зависимости от коэффициентов упругости связей геометрической модели. 
Общее решение задачи. Рассмотрим одноосное растяжение материала вдоль вертикальной оси 
силой Р.. В этом случае формулы (1) упрощаются, длины связей в текущей конфигурации: 


При одноосном растяжении сечение модели в перпендикулярной плоскости не искажает 
своей формы М№ = Во . Это дает возможность явно определить д, : 


1-9, = (а -9,)\2 
0 п (а ь)/2 = 5, =6,\2. (10) 


Кроме того, упрощается первая формула (4): 





216, +67 =246,+6.2-2а0, +6». (11) 
Углы в актуальной конфигурации: 
. В 40 — 0, ЙА 
зтф=—=————® и=-, 
Е 1+9, 4 
А + 
ера НОВ 9=ф, и 
Е 1-9 
40 — дь 
1 НЕЕ 
ыы А ато, 
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Уравнения статики (6) с учетом трансверсальной изотропии примут вид: 


26 зтф+ 2,6, +6, =0, (13) 
Е, =4К0, с0$ф + К ди. (14) 


Подставляя в выражение (13) формулы (10) и (12) после очевидных преобразований, по- 
лучаем: 


9.200: 6.=2 6: (15) 
где комплекс Со 20 (16) 
(2, +, + (2, +, +2А) 
Подставив соотношение (15) в (11), после преобразования приходим к уравнению относи- 
тельно д, с учетом зависимости (16): 


521-46 2а? СР (6,)) + 2а6/ (2 +2 С, (6,)) 246, -5.2=0. (17) 


Таким образом, при заданном 9, из нелинейных соотношений (16) и (17) определяется 





9, а далее по формулам (15) 0, и О Формулу для усилия (14) можно преобразовать таким 
образом: 
4-0, 





Е. =44б, 30: (18) 


ата 


4+9, 

Введем безразмерные величины: геометрические параметры приведем к начальной длине 

связи (4, а механические параметры - к коэффициенту упругости этой же связи Ко. Далее под 
старыми обозначениями будем иметь в виду новые безразмерные величины: 


1 д д д д [1 [4 [1 
И 0 Ее пб ПЕ Ве: СВЕ Е 
а Ч 40 90 40 а а а 
В такой постановке уравнения (15) - (18) принимают вид системы: 
2 


СЕЙ 
(о) М2, +) +9, + +24) 


521-С2)+26(/2+С)-26,-6/2=0 
бь=С,9; 
6, =6,\/2 


(19) 


Е. = | +4 5 | Кац. 
1 


Система (19) нелинейная, аналитическое решение ее не представляется возможным. Тем 
не менее, численное решение системы в зависимости от параметра удлинения д, производится 


без особых затруднений средствами систем компьютерной математики, например, пакета 
МаСаа. 

В результате получены нелинейные зависимости: 0,(0„), ди(д.), 6,(0.), Е,(9„). 
Частичная линеаризация задачи по деформациям (приближенное решение). Систему 


(19) можно упростить с помощью частичной линеаризации задачи по деформациям. Будем рас- 
сматривать малые деформации, когда можно пренебречь квадратами удлинений, но только в вы- 


ражении (11) с учетом безразмерных соотношений: 0, =0, —дь. 
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0. -—0 
Так как о ВЕСИ (20) 
0 [ +19 
Подставив (20) в (13) и (14), получим: 
2, 6 Е, +6, =0 
2+9, -0, (21) 
1+9 
Е, =(46 (6. -0,)— +0, Кац. 
а ( (9 5—5, в а@0 


В отличие от (18) в этом случае продольная деформация д,(0,) находится из квадрат- 
ного уравнения, а остальные характеристики удлинений и силы 0,(0.), д„(б„), Е,(б„) мож- 


но рассчитать аналитически с помощью несложных формул. 
На рис.4 графики представлены при А.=1 , Аь=1, К=2 и К„=1 для точного и приближен- 
ного решений. 
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Рис. 4. Сравнение деформационных и силовых характеристик модели 
для точного и приближенного решений 


На рисунке видно, что наибольшая погрешность между точным (нелинейным) решением и 
приближенным («полулинейным») наблюдается для величины [9 ‚ и достигает 30 % при дефор- 
мации порядка 1000 %, т.е. при десятикратном увеличении длины. При меньших деформациях 
порядка 10 % расхождение между нелинейным и линейным решением [9 значительно меньше -— 
около 3 %. Для других величин 0,(0,), 9,(6„), Е,(6,) расхождение даже при десятикратном 
удлинении не превышает 8 %. 
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Отметим, что параметр О, определяет характер влияния продольной деформации на по- 
перечные, то есть является параметром, отвечающим за «коэффициент Пуассона» модели. 
После определения зависимости силы от удлинения, поделив силу Ё` р На площадь попе- 


речного сечения до деформации 5 =44,”, получаем напряжение о’, а безразмерная 0, уже 
является относительной деформацией &. Далее можно вычислить нелинейную характеристику, 


О’ 
аналогичную линейному модулю Юнга Ё = — . Кроме того, из первого уравнения системы (21) 
[2 


можно найти нелинейную характеристику, аналогичную коэффициенту Пуассона у = 9%. 

а 
Полная линеаризация задачи (линейная задача). Получим коэффициент Пуассона и модуль 
Юнга в рамках линейной теории упругости при малых деформациях. Тогда д, —> 0, д, 0, 


и система (21) принимает вид: 
— р 6,) с (2^„ Е К,)0ь =0 


Е, =(20(6,-6,)+6,) Ка, . 
Далее приходим к соотношениям: 


[‹ 
6, => Г 6. 
К -+2К, + К, (22) 
2 
Е. = лы б.Каао . 
К +2, +К, 
Таким образом, коэффициент Пуассона трансверсально-изотропной модели в линейном 


гы К 
+2, + 
2 (ОК. +1 
Ен = лы Ка, /5 = (1+2. 
К -+2К, +А, 
Выбирая различные параметры модели А», А и К„ при нормировании А.=1, можно добить- 


ся различных значений коэффициента Пуассона и модуля Юнга. В табл.1 приведены значения 
коэффициента Пуассона и модуля Юнга, полученные при вариации значений параметров модели. 


случае: у ‚ а модуль Юнга в продольном направлении: 


К, +1))- 40/5. 


лин 
























































Таблица 1 
Влияние параметров модели на коэффициент Пуассона и модуль Юнга 

КЬ К Кл Улин Елин 

0 0 0 0 1 

1 1 1 0,250 0,625 
Е: р. 1 1 0,200 0,650 
то 3 1 1 0,167 0,667 
В 5 1 1 0,125 0,688 
$ 10 п 1 0,077 0,712 

20 1 1 0,043 0,728 
= 1 1,5 п 0,333 0,750 
Е 1 1 0,400 0,850 
= 1 2,5 1 0,455 0,932 
с 1 Е] 1 0,500 1,000 
в 1 5 1 0,625 1,188 
ы 1 1 2 0,167 0,667 
Е 1 1 3 0,125 0,688 
пх п 1 5 0,083 0,708 
2 1 1 10 0,045 0,727 
с 1 1 20 0,024 0,738 
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Данные значения позволяют понять роль каждого параметра модели. Как видно, коэффи- 
циент Пуассона зависит только от коэффициентов А, К, и к,. Увеличение упругости связи К 
влечет за собой увеличение коэффициента Пуассона. Увеличение же упругости связей А, и (К, 


влечет за собой уменьшение коэффициента Пуассона в обратно пропорциональной зависимости. 
И это понятно, ведь связи л и Ь препятствуют развитию деформации в перпендикулярной к на- 
правлению растяжения плоскости, уменьшая таким образом коэффициент Пуассона. 

При некоторых отношениях параметров модели можно получить значения коэффициента 


Пуассона, большие 0,5. Это возможно для анизотропного материала, когда К велико по сравне- 
нию с остальными параметрами А, и К,. 

Коэффициент Пуассона не зависит от величины коэффициентов А, модели, а только от их 
отношения. Модуль Юнга, наоборот, зависит как от отношения коэффициентов, так и от величи- 
ны этих коэффициентов. За размерность силы «отвечает» коэффициент АК а: Гаким образом, под- 


бирая параметры модели, можно получать разнообразные модули упругости и коэффициенты Пу- 
ассона практически существующих материалов. 

Отметим, что при значениях всех жесткостей пружин равных 1 коэффициент Пуассона ра- 
вен 0,25, что соответствует значению этого коэффициента, наиболее часто принимаемому в од- 
ноконстантной теории упругости. 

Уравнения модели для изотропного материала. В полностью изотропном материале 


ау = Ь, =с и Г = Ио = Ри. «Механическая» изотропия выражается в том, что 
К =А, = т И К, = К, = К. Тогда нетрудно показать, что система уравнений (19), описываю- 


щая одноосное растяжение вдоль связей а, примет вид: 


2 
С(брЕ 1 
22 +1) + 5,(4, +1) 
52(1-С2)+26,(/2+С,)-26,-6,2=0 





(23) 
бь=С1б1 
б,=6,\2 
Е ет 


2+ 6 
В линеаризованной форме при бесконечно малых деформациях получаем явное решение: 





6» = ы 0. =У:0, 

ЗК +1 (24) 

= И -бКаао = (1+2у(2%, +1))-6.К,ао. 

ЗК +1 
Здесь у = т является коэффициентом Пуассона. Модуль упругости Юнга 
а 
Е=| 1+ Ро К. а0/5 Е (1 +2у(2% +1))- 40/5 . Два размерных параметра Ка и 4 
1+ 


определяют размерность модуля упругости [А 1=НИМ, [а.]=м, [Е]=Па. В изотропном случае 
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Рассматриваемая модель представляет свойства известных линейно-упругих материалов. 
В табл.2 приведены значения параметров модели, соответствующие некоторым материалам. При 
этом выделяются две группы материалов: с коэффициентом Пуассона 0 <у<1/3 им>1/3З. Ин- 
тересно, что группа материалов с коэффициентами Пуассона у > 1/3 (сюда попадают, например, 
медь, серебро, свинец, золото) имеет повышенную пластичность, но небольшие модули Юнга [4]. 





















































Таблица 2 
Значения параметров модели для разных материалов 

Материал К) Ка, ГН/м у Е, ГПа 
Сталь 1,75 234 0,28 206 
Стекло кварцевое 0,347 19,0 0,17 7,5 
Стекло простое 1,00 112 0,25 70 
Кальций 4,43 1,46 0,31 2,6 
Цинк 1,42 130 0,27 100 
Железо 2,23 197 0,29 205 
Никель 8,00 70,7 0,32 210 
Хром 4,43 167 0,31 297 
Олово 33,0 3,98 0,33 45 
Вольфрам 3,00 269 0,30 350 
Бериллий 0,105 1005 0,08 300 
Серебро -7,00 -35,6 0,35 72 
Медь -3,36 -136 0,37 110 
Золото -1,62 -371 0,42 81 
Свинец -1,38 -119 0,44 16 

















При рассмотренном выше выборе значений жесткостей связей коэффициент Пуассона на- 
ходится в пределах 0 <у <1/3 (первая группа материалов). Для описания свойств материалов 
второй группы (у > 1/3) имеются две возможности: задать отрицательное значение жесткости А 
(безразмерное значение А,>-1 обеспечивает условие 0<у<1/2) или выбрать различные жест- 
К 


кости пружин связей при растяжении и сжатии. Во втором случае у“ и 
К +2, + К 

1 

модель будет иметь больше — свободных — параметров. При условии у<— 

2 


29 < +2 +А > Ю<2К, +А,. 
Для материалов с у >1/З3 (табл.2) связи модели предварительно деформированные: в 
размерном виде А, < 0, А, < 0, А, > 0, ав безразмерном А, > 0, А, <0. 


Надо заметить, что отрицательное значение жесткости пружин Ка и Кьне является фи- 


зически невозможным. Отрицательные значения К и А можно получить, например, следующим 


способом: к предварительно сжатым диагональным связям крепятся ненапряженные вертикаль- 
ные и горизонтальные связи. После освобождения всей конструкции от удерживающих сил диаго- 
нальные связи, стремясь вернуться в недеформированное состояние, разожмутся и растянут вер- 
тикальные и горизонтальные. В такой пространственной конструкции часть связей (диагональ- 
ные) будут находиться в сжатом состоянии, а оставшиеся связи (вертикальные и горизонтальные) 
будут растянуты. Таким образом, в полученной модели до определенного момента (до возврата в 
недеформированное состояние) горизонтальные связи не будут оказывать сопротивления попе- 
речному сжатию, а, наоборот, способствовать ему, что и будет соответствовать отрицательному 


коэффициенту жесткости пружин а и Б. 
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Заключение. В работе получены представления о свойствах математической модели нелинейно- 
упругих тел и о влиянии на них геометрических и механических параметров модели. Предложен- 
ные в статье определяющие соотношения, возможно, и сложные по виду формулы, достаточно 
естественны по смыслу параметров. Для разных частных задач путем последующих аппроксима- 
ций можно получить простые аналитические зависимости. 

При анализе свойств развиваемой математической модели основное внимание уделено 
двум нелинейным характеристикам, которым в линейной теории упругости соответствуют модуль 
упругости Юнга и коэффициент Пуассона. 


Библиографический список 

1. Лурье А.И. Нелинейная теория упругости / А.И. Лурье. - М.: Наука, 1970. - 940 с. 

2. Азаров А.Д. Моделирование больших деформаций полимеров / А.Д. Азаров // Совре- 
менные проблемы механики сплошной среды: тр. 1Х междунар. конф. - Ростов н/Д: ООО «ЦВВР», 
2005. - Т.1. - С. 13-16. 

3. Белл Дж. Ф. Экспериментальные основы механики деформируемых твердых тел. Ч.1. 
Малые деформации / Дж. Ф. Белл. - М.: Наука, 1984. - 600 с. 

4. Сандитов Д.С. Коэффициент Пуассона и пластичность стекол / Д.С. Сандитов, 
В.В. Мантатов, Б.Д. Сандитов // Журнал технической физики. - 2009. —Т. 79, выл. 4. - С. 150-152. 

5. Сандитов Д.С. Ангармонизм колебаний решетки и поперечная деформация кристалли- 
ческих и стеклообразных твердых тел / Д.С. Сандитов, В.В. Мантатов, Б.Д. Сандитов // Физика 
твердого тела. - 2009. — Т. 51, вып. 5. - С. 947-951. 

6. Фейнман Р. Фейнмановские лекции по физике. — Вып. 7. Физика сплошных сред 
/Р. Фейнман, Р. Лейтон, М. Сэндс. - М.: Мир, 1966. - 290 с. 


Материал поступил в редакцию 21.21.10. 


ВеГегепсе$ 

1. ге А.Т. Мейпетауа {еопуа иргидоз$Н / А.Т. Шиге. - М.: МацКа, 1970. -— 940 5. 
— п Кизчап. 

2. Агагоу А.О. Модейгоуате Бо’ деГогтасй ромтегоу / А.О. Ахагоу // Зоугетеппуе 
ргоМету тевап К! рМо$Нпо! эгеду: г. [Х теадипаг. КопЕ. — Козоу п/О: ООО «С\\В», 2005. - Т.1. 
— 5. 13-16. — м Кизчап. 

3. Вей р. Е. ЕКзреттета/пуе озпоуу теНап! аеГогтииетуй \\егауН *е|. СВ.1. Мауе 
ЧеГогтасй / О. Е. Ве. - М.: МацКа, 1984. - 600 $. — п Визчап. 

4. бапаКоу 0.5. Коейаепе Ричаззопа {1 р!азисппо$Ё $еКо| / 0.5. бапайом, У.\У. Мапаюом, 
В.О. бапаКо\ // Зигпа! {ерпспезКо! Я71. - 2009. - Т. 79, мур. 4. - $. 150-152. — Тп Кизчап. 

5. бапаКом 0.5. Апдагтопят КоеБапй гезпек! 1 рорегесппауа ЧЕРогтадуа Ки$аШспезкИ | 
$еКюоБгагпуй \егауй *=!| / 0.5. бапаКом, У.\. Мапкаюу, В.О. бапаКо\ // АжКа \уегдодо {а1а. - 2009. 
—Т. 51, мур. 5. - 5. 947-951. — т Визчап. 

6. Геттап В. Реттапоу$Ке 1еКсИ ро Я7Же. - Мур. 7. Е2Жа зрюзйпуН эге4 / К. Геттап, 
В. ГеКоп, М. 5епд$. - М.: Ми, 1966. - 290 $. — Тп Кизчап. 
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(Кееагсн пище оЁ Меспапс$ апа АррИе Маетачс$ аКег 1.1. Могомсп), 
О.А. АХАКО\У 

(Роп За ТесНтса! Упмег®Ку) 


А 3 тесватса! тодЕ! Гог аезспирноп ог 1агде дегогтаНоп$ о! поп-Шпеаг е/азИс Бобез [5 обегец. Тре тодЕ!/ 
Че5спе5 уапоиз {урез ог ап5оору, 5исй аз оппогор, гапзуегсе/ /5о{торг, аз ие! а5 15оторс Бобе5. Тре 
сазе оЁ итажа!/ епзюп 15 сопз/аегед. Тиго спагасепйуНс$ ог {те поппеаг таепа! хтйаг юЮ {те [пеаг Уоипда’5 
тоаШ/из апа Ро/55оп'5 гаНо аге туезИдае4. Тре Гппег-тодЕ! рагатеегз ейесйпу таепа/ зресйсайоп5 аге 
апа/угец. 

Кеуигога5: поп/теаг е/а5Ёсйу, аг/зо{гору, та етайса! тодЕ!/ Уоипд’5 тоди/и<, Рой55оп‘5 гайд, итажта! {еп5юп. 
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